粘性流体大幅晃动的ALE有限元模拟

采用有限元法数值求解了具有自由液面大幅移动边界的Navier－Stokes方程。对流体区域采用了任意的拉洛朗一欧拉（ALE）运动学描述，网格结点可以任意移动而不依赖于流体的运动，结合了拉格朗目描述易于处理移动边界与区域变形的优点和欧拉描述可克服单元缠结的优点，提出了简单合理的网格更新方法，能精确地跟踪运动的自由液面。为了更精确地处理强对流项的作用，采用了迎风流线Petrov－Galerkin（SUPG）加权余量法建立有限元方程。算例表明本方法对贮箱内流体的非稳态大幅晃动过程的数值模拟非常成功。     

二维非结构网格Euler方程的LU-SGS算法

研究了将二维非结构三角形网格重新排序后，用LU-SGS隐式算法计算Euler方程的方法。分析了网格不平衡情况下该方法的可行性和应用效果。计算了不同翼型在不同流动条件下的二维Euler方程，并与四步Runge-Kutta显式方法进行了比较，结果验证了此方法具有良好的计算和收敛效果，以及在粘性计算方面的潜力。运用此方法对某预研型号的头部外形进行了设计和选择。该方法克服了以往隐式方法大量耗费内存的弱点，达到了计算耗时短和占用内存少的统一。     

圆管内不可压粘流的分布式并行算法

用并行算法求解了圆柱通道内不可压层流流动。采用有重叠的条形区域分裂法时计算域进行分区，分别求解动量方程和压力校正方程，对压力校正方程采用块修正法来加速线代数方程的收敛速度。计算结果表明有较好的并行加速比和并行效率。     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

A coupling VWM/CFD/CSD method for rotor airload prediction

A coupling fluid-structure method with a combination of viscous wake model(VWM),computational fluid dynamics(CFD) and comprehensive structural dynamics(CSD) modules is developed in this paper for rotor unsteady airload prediction. The hybrid VWM/CFD solver is employed to model the nonlinear aerodynamic phenomena and complicated rotor wake dynamics;the moderate deflection beam theory is implemented to predict the blade structural deformation; the loose coupling strategy based on the ‘delt method' is used to couple the fluid and structure solvers.Several cases of Helishape 7A rotor are performed first to investigate the effect of elastic deformation on airloads. Then, two challenging forward flight conditions of UH-60 A helicopter rotor are investigated, and the simulated results of wake geometry, chordwise pressure distribution and sectional normal force show excellent agreement with available test data; a comparison with traditional CFD/CSD method is also presented to illustrate the efficiency of the developed method.     

用粘性涡方法计算伞状物绕流

文章用粘性涡方法计算了二维伞状物绕流问题。首先阐述了涡方法的基本理论 ,包括涡元的对流、粘性扩散、物面涡元的脱落、物体所受气动力计算等 ,然后介绍了使用涡方法的一般步骤。最后将完全充满的伞衣简化为二维圆弧壳 ,计算了伞衣内外压力分布、阻力系数及升力系数     

Local stability analysis of interface region of astrophysical viscous shear flows with a gradual velocity gradient

In this paper, we extend the stability analysis of cold sharp shear flows to warm astrophysical cases with, inevitable, gradual velocity gradient in the interface region in the presence of viscosity effect. Using linear perturbation theory as well as the local approximation method, the instability growth rate of the excited electromagnetic modes has been investigated for the relativistic and non-relativistic cases of solar wind interacting with interstellar plasma medium. Results show that astrophysical shear systems with a small velocity gradient in the transition region are more stable rather than larger ones. Moreover, dependent on the viscosity coefficient value, the viscosity effects could have a positive role on the instability growth rate of the system in some range of initial bulk velocity, while it plays a destructive role in other velocity ranges.     

高超声速乘波飞行器气动特性研究

用计算流体力学和风洞试验的方法对以锥导乘波体为基础生成的高超声速乘波飞行器的气动性能进行了研究。结果表明，以马赫数6，攻角4度为设计状态的乘波体，在马赫数5～7，攻角4～6度的范围内，都具有良好的气动特性，升阻比接近4。最后，提出了一个简单的以参考温度方法为基础的粘性阻力分析方法。该方法配合使用风洞试验和计算流体的结果，可以用来验证计算流体中难以计算准确的粘性阻力，也可以用来分析在风洞试验难以直接得到的粘性阻力。对于工程上的粘性阻力分析是一个有用的办法。     

Development of a coupled supersonic inlet-fan Navier–Stokes simulation method

A coupled supersonic inlet-fan Navier–Stokes simulation method was developed by using COMSOL-CFD code. The flow turning, pressure rise and loss effects across blade rows of the fan and the inlet-fan interactions were taken into account as source terms of the governing equations without a blade geometry by a body force model. In this model, viscous effects in blade passages can also be calculated directly, which include the exchange of momentum between fluids and detailed viscous flow close to walls. NASA Rotor 37 compressor test rig was used to validate the ability of the body force model to estimate the real performance of blade rows. Calculated pressure ratio characteristics and the distribution of the total pressure, total temperature, and swirl angle in the span direction agreed well with experimental and numerical data. It is shown that the body force model is a promising approach for predicting the flow field of the turbomachinery. Then, coupled axisymmetric mixed compression supersonic inlet-fan simulations were conducted at Mach number 2.8 operating conditions. The analysis includes coupled steady-state performance, and effects of the fan on the inlet. The results indicate that the coupled simulation method is capable of simulating behavior of the supersonic inlet-fan system.     

Mach reflection in steady supersonic flow considering wedge boundary-layer correction

Mach reflection in steady supersonic flow is an important phenomenon having received extensive studies, among which simplified theoretical models to predict the size of Mach stem and other flow structure are of particular interest. Past efforts for such models were based on inviscid assumption while in real cases the flow is viscous. Here in this paper we consider the influence of wedge boundary layer on the Mach stem height. This is done by including a simplified boundary layer model into a recently published inviscid model. In this viscous model, the wedge angle and the trailing edge height, which control the Mach stem height, are replaced by their equivalent ones accounting for the displacement effect of the wedge boundary layer, with the boundary layer assumed to be laminar or fully turbulent. This viscous model is shown to compare well with numerical results by computational fluid dynamics and gives a Mach stem height as function of the Reynolds number and Mach number. It is shown that due to the viscous effect, the Mach stem height is increased, through increasing the effective wedge angle.